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Abstract 

We have investigated the kagome ice state in the frustrated pyrochlore oxide Dj^T^Oy 
under magnetic field along a [111] axis. Spin correlations have been measured by 
neutron scattering and analyzed by Monte-Carlo simulation. The kagome ice state, 
which has a non-vanishing residual entropy well established for the nearest-neighbor 
spin ice model by the exact solution, has been proved to be realized in the dipolar- 
interacting spin ice Dy2Ti207 by observing the characteristic spin correlations. The 
simulation shows that the long-range interaction gives rise to only weak lifting of the 
ground state degeneracy and that the system freezes within the degenerate kagome 
ice manifold . 
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The spin ice model, an Ising ferromagnet with a local <111> easy axis on the 
pyrochlore lattice interacting with nearest-neighbor exchange coupling, has 
attracted much attention because of its intriguing frustration mechanism pQ. 
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It has macroscopically degenerate ground states characterized by the "two-in 
and two-out" local structure on each tetrahedron, equivalent to the "ice rule" 
in the water ice [2]. The spin ice behavior, e.g. the residual entropy [3J, was 
discovered in the rare earth pyrochlore oxides, Ho 2 Ti 2 7 [I], Dy 2 Ti 2 7 [3] 
and Ho 2 Sn 2 07 j5]. However in these systems Ising spins on the rare earth sites 
interact via a long-range dipolar interaction, and the mechanism showing the 
spin ice behavior in these dipolar spin ice was puzzling. A numerical account 
was provided by the Monte-Carlo (MC) simulation studies pQ, showing that the 
ground state manifold is approximately preserved in the dipolar model. More 
recently an elegant explanation using an analytical method was proposed in 
a theoretical work [TJ. 

The macroscopic degeneracy of the spin ice model is partly or fully lifted by 
magnetic fields [8j. Along a [111] axis the pyrochlore lattice can be viewed as 
an alternating stacking of kagome and triangular layers. For field along the 
[111] axis, the degeneracy is partly lifted, where random spin configurations are 
retained only in the kagome planes, and has been studied by an exactly solvable 
dimer model [91TT0] . Surprisingly, this kagome ice state based on the nearest- 
neighbor interacting model, was suggested to be observed experimentally in 
the dipolar spin ice compound [TT][T2"] . 

Unlike the case in zero field, the effects of the long-range dipolar interaction 
have not been clarified in magnetic field, and hence it is still puzzling whether 
the degeneracy of the kagome ice manifold is preserved in the dipolar spin 
ice. To resolve this question, we have measured static spin correlations of 
the dipolar spin ice compound Dy 2 Ti 2 7 in magnetic field along the [1 1 1] 
axis by neutron scattering and observed characteristic spin correlations. By 
comparing with MC-simulation results based on the dipolar model, we have 
demonstrated that the kagome ice state really occurs in the dipolar spin ice. 

Neutron scattering experiments were performed using the triple-axis spec- 
trometer GPTAS installed at JRR-3M JAERI. Two single crystalline samples 
of Dy 2 Ti 2 7 , which were grown by the floating-zone method using an infrared 
furnace, were used for the experiments. The samples were mounted in a di- 
lution refrigerator. MC-simulations were carried out using the facilities of the 
Supercomputer Center at ISSP University of Tokyo, which utilize the standard 
Metropolis single spin flip algorithm. 

We show neutron intensity maps in the scattering plane perpendicular to the 
[1 1 1] axis in zero and finite(if = 0.5 T) fields at T = 0.3 K obtained from the 
experiments and MC-simulations based on the dipolar model in Fig. [TJ Ex- 
perimental results show definite differences between the spin correlations at 
H = 0.0 and 0.5 T, as shown in Fig. [1] (I-a) and (II-a) respectively. The char- 
acteristic correlations at Q = (2/3, —2/3, 0)(Q - corre l a tion), the positions 
are represented by black circles in (II-a), is found at H = 0.5 T. 
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Q = (2/3 -2/3 0) h in (h -h 0) (r.l.u.) 
Intensity n , ■ 
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Fig. 1. Neutron intensity maps in the perpendicular plane to the [1 1 1] axis. Ex- 
perimental results (I-a) at T = 0.3 K and H = 0.0 T and (Il-a) at T = 0.3 K and 
H = 0.5 T. MC-simulation results (I-b) at T = 0.3 K and H = 0.0 T and (Il-b) 
at T = 0.3 K and H = 0.5 T. Energy parameters used for the simulations in the 
dipolar model are J x = -3.72 K, J 2 = 0.1 K, J 3 = -0.03 K and D = 1.41 K, where 
the 1st, 2nd, 3rd nearest neighbor exchange interactions and the dipolar coupling 
constant respectively. The values of J%, J3 and D are assigned to those reported 
previously |6|13j and the value of J2 is determined to reproduce the temperature 
dependences of the specific heat in finite fields [T4"]. 

The MC-simulations based on the dipolar model well reproduce the neutron 
intensity map at H = 0.5 T(Fig. Q] (II-b)), as well as that at H = 0.0 T(Fig. CD 
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(I-b)). The characteristic Qq- correlation observed in the experiment is clearly 
reproduced in Fig.[TJ(II-b). The (3o- corre l a ti°n was proposed by the theoretical 
work based on the nearest-neighbor model [H], and hence, it is a good signature 
of the kagome ice state. The agreement between the experimental and the MC- 
simulation results on the static spin correlations, especially the observation of 
the Qo-correlations, indicates that the degeneracy of the kagome ice manifold 
is approximately preserved in the dipolar spin ice Dy2Ti2C>7. Spin correlations 
at H = 0.5 T do not show significant temperature dependence below 0.75 K 
in both the experiments and the single spin flip MC-simulation. It strongly 
suggests that the system freezes within the degenerate kagome ice manifold. 

In summary, we have presented the investigation of the neutron scattering 
experiments of the dipolar spin ice compound Dy2Ti2C>7 in the magnetic field 
along the [1 1 1] axis. The analysis of the experimental results by the single 
spin flip MC-simulations reveals that the kagome ice state is realized in the 
dipolar spin ice. 



References 



[I] S. T. Bramwell and M. J. P. Gingras, Science 294 (2001) 1495. 
[2] P. W. Anderson, Phys. Rev. 102 (1956) 1008. 

[3] A. P. Ramirez, A. Hayashi, R. J. Cava, R. Siddharthan, B. S. Shastry, 
Nature (London) 399 (1999) 333. 

[4] M. J. Harris, S. T. Bramwell, D. F. McMorrow, T. Zeiske, K. W. Godfrey, 
Phys. Rev. Lett. 79 (1997) 2554. 

[5] H. Kadowaki et al., Phys. Rev. B 65, 144421 (2002). 

[6] B. C. den Hertog and M. J. P. Gingras, Phys. Rev. Lett. 84 (2000) 3430. 

[7] S. V. Isakov, R. Moessner, S. L. Sondhi, Phys. Rev. Lett. 95 (2005) 217201. 

[8] M. J. Harris, S. T. Bramwell, P. C. W. Holdsworth, J. D. M. Champion, Phys. 
Rev. Lett. 81 (1998) 4496. 

[9] R. Moessner and S. L. Sondhi, Phys. Rev. B 68 (2003) 064411. 

[10] M. Udagawa, M. Ogata, Z. Hiroi, J. Phys. Soc. Jpn. 71 (2002) 2365. 

[II] K. Matsuhira, Z. Hiroi, T. Tayama, S. Takagi, T. Sakakibara, J. Phys.: Condens. 
Matter 14 (2002) L559. 

[12] Z. Hiroi, K. Matsuhira, S. Takagi, T. Tayama, T. Sakakibara, J. Phys. Soc. Jpn. 
72 (2003) 411. 

[13] J. P. C. Ruff, R. G. Melko, M. J. P. Gingras, Phys. Rev. Lett. 95 (2005) 097202. 
[14] Y. Tabata et al., |cond-mat/0607808. 



4 



